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Abstract 

We consider a scenario where our universe is taken as a three-dimensional domain 
wall embedded in a five-dimensional Minkowski space-time, as originally proposed by 
Brito and collaborators pQ . We explore the existence of a richer class of soliton solutions 
and their consequences for the acceleration of the universe driven by collisions of bulk 
' particle excitations with the walls. In particular it is shown that some of these solutions 

£T) ■ should play a fundamental role throughout the expansion process. 
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1 Introduction 



Very important systems described by quantum field theories are intrinsically nonlinear; the 
Standard Model and the Quantum Chromodynamics are classical examples. By exploring 
deeply such systems or effective nonlinear models, even at the classical level, one has shown 
the increasing importance of the soliton solutions (classical solutions with finite and localized 
energy) and their broad applications [2j-[5j. Soliton solutions in quantum field theory de- 
scribe, for instance, monopoles, magnetic vortices, instantons in quantum chromodynamics, 
cosmic strings and magnetic domain walls. Finding exact classical solutions, particularly 
solitons, is one of the problems on nonlinear models with interacting fields. When one has in 
hands a systematic method, as the one offered by Rajaraman [BJ, this task becomes easier, 
even though one has to explore the consequences of the classical solutions, as well as their 
possible realizations in the nature. As pointed out by R. Rajaraman and E. Weinberg [TJ, 
in such nonlinear models more than one time-independent classical solution can exist and 
each one of them corresponds to a different family of quantum states, which come into play 
when one performs a perturbation around those classical solutions. 

The method in [6_], usually called trial orbits method, is a very powerful one presented for 
finding exact soliton solutions for nonlinear second-order differential equations of models 
with two interacting relativistic scalar fields in 1+1 dimensions, and it is model indepen- 
dent. The trial orbits method have been applied to the special cases whose soliton solutions 
of the nonlinear second-order differential equations are equivalent to the soliton solutions 
of first-order nonlinear coupled differential equations, the so-called Bolgomol'nyi-Prasad- 
Sommerfeld (BPS) topological soliton solutions [8]. Some years ago one of us presented a 
method for finding additional soliton solutions for those special cases whose soliton solu- 
tions are the BPS ones [9] . Latter, that approach was extended by considering more general 
models [TU]. Furthermore, that method shows how to get the general equation of the orbits, 
explains how the different solutions connect the different vacua of the model under analysis 
and, as a novelty, presents a class of soliton solutions with a kink-like profile for both fields 
with its minimum energy (BPS energy) smaller than that of the usual solution which exhibits 
a kink-like configuration for one of the fields and a lump-like configuration for the other one. 
Moreover, the stability of the quantum states corresponding to these new soliton solutions 
can be shown on the same basis presented in [TTj . 

The BPS soliton solutions have found applications in a great variety of natural systems 
whose dynamics can be approximately described by nonlinear quantum field models for 
interacting scalar fields in 1+1 dimensions [12]. Those kind of models have been generalized 
by including into the Lagrangian density some minimal terms that break Lorentz and CPT 
symmetries [13J. By analyzing the consequences of additional soliton solutions in a given 
nonlinear model, some of us have shown [T3] , by using the method developed in [5], that those 
nonlinear Lorentz breaking models in 1+1 dimensions exhibit additional soliton solutions 
whose BPS energies are smaller than those found in [13], and that more general Lorentz 
breaking models in 1+1 dimensions, which admits soliton solutions, can be built. 

In this work, we explore more deeply the classical solutions found in [9], particularly a 
class of degenerate soliton solutions [T5| [TS] , for the nonlinear model of two interacting scalar 
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fields in 1+1 dimensions [12]. In special, we analyze the consequences that those additional 
soliton solutions bring for the scenario of accelerating universes. This scenario was recently 
conceived by Brito and collaborators [T] and it is within the context of the extra-dimensions 
[TB]- [24J . Our analysis is done by following a similar approach to that of the reference [T|. 
It is shown that, for a critical value of the degeneracy parameter, the reflection coefficient 
of the bulk particles over the wall is maximized. That implies into a degree of expansion of 
the universe which is compatible with the present observed acceleration. 

In the second section we present the model we are going to work with and discuss the 
influence of each set of soliton solutions, in the context of accelerating universes driven by 
the quantum states (bulk particles). In the third section we present final comments about 
the impact of those new solutions over the expanding process. 

2 Domain walls and accelerating universes 

A very important and intriguing modern physical problem is that of finding a way to explain 
the observed accelerated expansion of the universe. On the other hand, the recent cosmo- 
logical data indicates that a relevant part of the energy of the universe would be a kind of 
dark energy [25]. In fact, that dark energy is supposed to be the responsible for that accel- 
eration. As a consequence, many authors look for a deep understanding of these subjects. 
A very interesting possibility is the one associated to the so-called brane worlds [24] . In a 
recent work, Brito and collaborators [1] conceived a scenario where our universe is taken 
as a three-dimensional domain wall embedded in a five-dimensional Minkowsky space-time, 
where the elastic collisions of bulk particles with the walls would be the ultimate reason for 
the universe acceleration. 

The model under analysis here is precisely the same one considered in the reference [I], 
which is the scalar sector of a five-dimensional supergravity theory obtained by means of 
dimensional compactification of a higher dimensional supergravity, where the Lagrangian 
density [TJ looks like 

C=\detg^!^~M 3 R 5 + G AB 

where 1/M is the five- dimensional Planck length, R§ is the Ricci scalar, and Gab is the metric 
on the scalar target space. Now, considering the limiting situation where M » 1 << Mpi, 
so that the five-dimensional gravity decouples from the scalars, one is left with a 3-brane 
which is a classical 3d domain wall embedded in a 5d Minkowski space [I]. Thus, in this 
limit, and choosing the case with two scalar fields in the target space, we are working with 
a nonlinear model with interacting scalar fields in 1+1 dimensions whose time-independent 
equations of motion (Eqs. below), are those one-dimensional static equations for the 
scalar fields taken into account in the reference [T]. We follow the same route as in that 
reference, but here we consider more general set of static solutions and compare our results 
with those obtained in [T] where the bulk particles, which are the quantum excitations around 
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the classical solutions, collide with a 3d domain wall described by a kink. The essential idea 
is to show that the situation is richer than analyzed in [I] , and that from a complete solution 
as the one we discuss here, important consequences for the expanding scenario show up. 

The model was introduced before to study BPS soliton solutions [?] and, recently, it was 
explored in the context of brane worlds [16]. It consists of two interacting real scalar fields 
in 1+1 dimensions and it is described by the Lagrangian density 

C^d^f + \{d,xf ~V{faxl (2) 
where the potential is given by 

V(fa x) = ^A 2 (0 2 - a 2 ) 2 + (2/i 2 + A ^)0Y - A /i a 2 x 2 + \^\\ (3) 

The distinctive property of this model is that its potential can be written in terms of a 
so-called superpotential as 

. i ( dw(4,, x ) \ 2 , i ( dwjfad Y m 



where the superpotential is 



2 i 2 

a ]+nx 



(5) 



Hence, finding the classical solutions with minimum energy for the time-independent equa- 
tions of motion 

fl = dV ^ fx = ^Y_ (6 ) 

dx 2 d(f) dx 2 dx ' 

is equivalent to find the classical solutions with minimum energy for the time-independent 
first-order differential equations 

0' = ^(0, X ), x' = W x (<f>, X ). (7) 

In the above equations the prime means the derivative with respect to the space coordinate, 
and W(f, (W x ) stands for the partial derivative of W(<f>, x) with respect to the 4> (%) field. 
The minimum energy (BPS energy) [8] for nonlinear systems described by the Lagrangian 
density (J2J) with potentials written as (j3J) is found to be given by 

E B ps=\W((t> j ,x j )-W(<i> i ,Xi) \ , (8) 

where fa and x% mean the i — th vacuum states of the model. 

Applying the method developed in [9], we note that it is possible to write the relation 
d(p/W^ = dx = dx/W x , where the differential element dx is a kind of invariant. Thus, one 
is lead to 
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d x w x w 

This is in general a nonlinear differential equation relating the scalar fields of the model. 
If one is able to solve it completely for a given model, the function cf) (x) can be used to 
eliminate one of the fields, rendering the equations (JZJ) uncoupled and equivalent to a single 
one. Finally, the resulting uncoupled first-order nonlinear equation can be solved in general, 
even if numerically. By substituting the derivatives of the superpotential (JSJ) with respect to 
the fields in (J9j), it can be rewritten as a linear differential equation, 

dp X , . 

/ = X, 10 

by the redefinition of the field, p = (f) 2 — a 2 . Now, the general solutions are easily obtained 

as 

p( X ) = <p 2 -a 2 = c X X/fl ~ X\ for A ^ 2p, (11) 

where Cq is an arbitrary integration constant and, for the sake of simplicity, we will restrict our 
study to the case where A 7^ 2p. Then, we substitute the above solutions in the differential 
equation (J7|) and obtain the following first-order differential equation for the field x{ x )i 
obtaining 

^ = ±2p X] J^ + c x x/ » ~ J-^X 2 > (12) 

Despite the fact that in general explicit solutions for the above equation can not be obtained, 
one can verify numerically that the solutions belong to the same classes, and some of those 
classes of solutions can be written in closed explicit forms. In those last cases we are able to 
obtain the several types of soliton solutions, whose consequences we discuss below. 

The set of solutions used by Brito et al [T], called here type-A kink [§] to distinguish it 
from other types of kink solutions we present for this model, can be obtained by means of 
the method described in the previous section just by taking cq = in the expression ffl~2]) . 
In this case that equation can be solved analytically for any value of A and p, in the range 
A > 2 p. We have the following solutions for x{ x ) an d <t> (x)'- 



Xa(x) = a\ — — sech(2yuax), 4>a(x) = ±a tanh(2/iax). (13) 

V P 

In this case the BPS energy is given by Eg PS = |Aa 3 . 

Other soliton solutions can be found when one considers the integration constant Co 7^ 
[15] . It was found in [9j [101 [HI [16] that in three particular cases the equation (fT2l) can be 
solved analytically. For Co < —2a and A = p it was found that the solutions for the x( x ) 
field are lump-like solutions, which vanishes when x — > ±00. On its turn, the field <fi{x) 
exhibits a kink-like profile. They also connect the vacua of the model and also have BPS 
energy E^ PS = |Aa 3 . These classical solutions can be written as 



v 2o - (1) ^cg-4sinh(2/iax) 

Xa( x ) = rr, and ( x ) = A » ( 14 ) 

— 4 cosh(2/iax) — cq Jcq — 4 cosh(2pax) — Co 
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An interesting aspect of these solutions is that, for some values of c < —2a, $^(x) 
exhibits a double kink profile. We can speak of a formation of a double wall structure. In 
Figure 1 we plot some typical profiles of the soliton solutions in the case where A = fi, both 
when Co is close to its critical value (co = —2 in this case) and far from it. Both fields are 
there represented. One can verify that the distance from one wall to the other one increases 
as Co approaches its critical value. For the critical value of Co the double wall structure 
merges into a single one. 

Similar behavior is also noted in the classical solutions for A = 4/i and Co < 1/16. In this 
case the fields profile look like 

&\x) = 2a — , (15) 



'y/l — 16co cosh(4/iax) + 1 
and 

ZWf \ ft TcT- sinh(4/iax) 

<t>A '{x) = VI - lQc a = = — — — . (16) 

y/l — 16c cosn(4//aa;) + 1 

In the above solution one can also see the double kink profile for some values of cq, and the 
increasing of the distance from one wall to the another as Co approaches its critical value 
(c = 1/16 in this case). Once again, at the critical value of c , the double wall structure 
coalesces into a single wall. In Figure 2 it is shown the behavior of the energy density 
for A = 4/i. From that, it is quite evident the appearance of the double walls when one 
approaches Co = l/16a 2 . 

Finally, very interesting analytical soliton solutions were shown to exist when one takes 
A = jj, and the critical parameter Co = —2a and for A = 4/i and the critical parameter 
c = l/16a 2 , in the equation f[T2"j) . The novelty in these cases is the fact that both, the 
x{x) field and the <ft(x) field present a kink-like profile and the BPS energy is half of that of 
type-A kinks, that is E§ PS = |Aa 3 . We call this set of solutions as type-B kinks. For A = \i 
and Co = —2a the classical solution for the fields can be shown to be given by 

Xb\ x ) = =t tanh (/ia a;)) and 4>b\x) = ^(tanh(/xaa;) =F 1)- (17) 

For c = l/16a 2 and A = 4/i, the following set of type-B kinks is obtained 

jSHx) = _^^if^)^(f^) and g) (x) = "(i Tt anh(2^)). (18) 

Jcosh(2/iax) ^ 

The type-B solutions also connect the vacua of the model, but this time one could interpret 
these solutions as two kinds of torsion in a chain, represented through an orthogonal set of 
coordinates (f) and x- m the plane (</>,x); the type A kinks correspond to a complete torsion 
going from (—a, 0) to (a, 0) whereas the type B corresponds to a half torsion, where the 
system goes from (—a, 0) to (0, a), in the case where (A = /x) for instance. 

Now, by proceeding as in PQ, we perform a linear perturbation of the x{ r i t) field around 
the classical solutions, that is 

X(r, t) = x(r) + C(r, t) and 0(r, t) = 0(r), (19) 
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where x(r) and 4>(r) are the classical solutions (background fields) and £(r, t) is the quantum 
field. By expanding the action up to quadratic terms in the quantum fields we obtain second- 
order differential equations for the quantum fields 

d fl d»( + V xx (r)C=0, (20) 

that is, the quantum field obeys a Klein-Gordon equation with an effective potential V xx (r) 
which is obtained by taking the second derivative of the potential given in (j3J) with respect 
to x( r ) an d evaluated at the classical solutions x(r) and 4>(r) 

V xx (x) = 2(2/i 2 + A /i)0 2 + 6/iV - 2A /i a 2 . (21) 

If we consider our model as a five- dimensional one and that the fluctuations can be 
expanded in terms of plane-waves in the space-time coordinates, that is 

C(r, t) = C(r) exp[-i(ut - k x x - k y y - k z z), (22) 

with £(r) a function of the fifth coordinate, we obtain the following Schrodinger-like equation 

(-^+^wjc(r)=eC(r), (23) 

where e = u 2 + k 2 x + k 2 + k 2 and we are considering that the static solutions described in the 
previous section are now functions of the coordinate r. 

Now, we compare the effect of all the above soliton solutions over the acceleration scenario 
proposed in pQ. For this, we begin by considering the static solutions given by the equations 
(JHj). By substituting them in the effective potential expressed in ( |2~T]) . we find that the 
Schrodinger-like equation for the quantum excitations is 



y dr 

where 



- 2 +V e f f - 1) (r))ar)=e((r), (24) 



V}ff 1] {r) = 6m V >' 4 ) sinh (2^r) + 4^ _ ^ (2g) 



' Cq — 4 cosh(2 fiar) — c 
and, when c = —2, the effective potential, for fi = A, becomes 



Vjff 1} (r) = n 2 a 2 (l + 3 tanh 2 (/xar)). (26) 



eff 

Note that, in the case analyzed in p], the effective potential is 

A 



V${r) = 4fia 2 -4fia 2 4- - ) sech 2 (2 fiar) . (27) 



In Figure 3 it is shown the behavior of the effective potential coming from the domain wall 
solutions in four situations: the case studied by Brito et al. and for the case of degenerate 
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solitons when Co is far from its critical value, near it and at the critical point (cq = —2). From 
that Figure one can perceive that, apart from the appearance of the two wells potential, for 
Co close to the critical value, the case where the reflection coefficient is bigger is probably 
the case studied in [TJ. 

The more remarkable result comes to life for A = 4/i. By substituting the equations (j!5p 
and (fTBl) in (1211) , one can verify that the quantum excitations of the x{ r i t) field satisfies the 
effective Schrodinger equation fl23|) with an effective potential given by 



T/ (A-2)/ v 1Q 2 2 (1 - 16co) sinh (4/xar) + 2 2 2 

K// V) = 12/x a — — -y - 8/i a , (28) 

I yl — loco cosn(4yU a rj + 1 J 



and, when Co = 1/16, the effective potential becomes 

V}ff 2 \r) = fi 2 a 2 sech 2 {2 fi a r) [2 + 5 cosh(4/iar) + 3 sinh(2 /x or)] . (29) 

It can be observed from equation [27] that, for A = 4/i, one has a constant effective potential 
V}ff (r) = Ana 2 , so there is no reflection of the bulk particles. This is an artifact of the 
particular solution considered in pQ. In contrast with this, the solutions presented here lead 
to effective potentials such that the reflection coefficient increases as Cq approaches its critical 
value. Those effective potentials are presented in Figure 4. One can note that, for cq far 
from and close to its critical value, there is a non-zero probability for the bulk particles to be 
transmitted, whereas, for the critical case, the effective potential acquires a ramp-like profile, 
so granting that the bulk particles with energy below the height of the ramp are certainly 
reflected by that "infinitely large barrier" . 

For a given amount of energy available to build up structures in the original universe, we 
argue that it should be distributed, for instance, like 

M 

£ t = J2 n{ ^bps, (30) 
i=i 

where stands for the number of defect structures of the i-th specimen and 

(') (') ^ 2 

£-bps = Ebps + 7r~5 (31) 
2m 

with E^ PS being the BPS energy of a given soliton configuration, and k r the momentum of 
the excitation with mass m hitting the wall. 

In order to estimate the force acting on the wall by the incoming excitations, we consider 
quasi-elastic collisions to guarantee the energy conservation in the process. In other words, 
the kinetic energy acquired by the wall comes from the hitting particles. In that case, the 
momentum transferred to the wall is given by 

AP N reflected {2 kr) = N TI R(2k r ), (32) 

where R is the reflection coefficient and Nti is the total number of incoming particles. From 
this, we can conclude that the acceleration of the wall is 
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wan , , A . 5 

(33) 

M too h At 



with 

2 k r N TI 
M wall At ' 

and is the flux of incident particles. 

From now on we follow the brane inflation scenario presented in pQ, and similarly in 
[21]. In that scenario one supposes that we live in a three-dimensional domain wall which 
is considered as the surface of a hypersphere, whose radius increases due to the hitting 
bulk particles. The increasing of the radius can be modeled as r(t) = a(t) f , where a(t) is 
the known scale factor and f is the constant radius of the bubble in the comoving system. 
Consequently, in this cosmological model, we have the metric of the 4d universe given by 

ds\ = dt 2 -a(t) 2 (dx 2 + dy 2 + dz 2 ). (34) 

From this point of view, the acceleration of the wall is proportional to the second-derivative 
in time of the scale factor, that is a wa u = d(t)r, and we have the following relation 

ii{t) = KR. (35) 

where k = a/r. 

It can be seen from Figure 4 that the effective potential barrier width increases when the 
degeneracy parameter Co approaches its critical value Cq. As a consequence, the reflection 
coefficient becomes asymptotically constant. Conversely, for cq far from the critical value 
the reflection coefficient tends to zero, since the width and height of the barrier diminish in 
this limiting situation. Thus, when c is close to the critical value, one has 

a(t)~^kt 2 , (36) 

which is compatible with the present observed behavior of the universe. Furthermore, in 
the opposite limit, when Co is far from the critical value, the reflection coefficient approaches 
zero. In this situation one can suppose that this happens exponentially, so the scale factor 
is given by 

as requested by an inflationary era. 

At this point it is important to remark that Cq is a function of the parameters /j, and A. 
This can be noted from the fact that its value changes when the relation among A and \i 
changes (see the cases analyzed above). This property allows us to suppose that Cq would 
have a time dependence, as well as the parameters of the potential. That dependence might 
have its origin in the universe expansion itself. The temperature of the universe decreases 
along with its expansion and, consequently, the parameters of the potential would exhibit a 
temperature dependence. Therefore, Cq approaches the value of cq the universe was created 
with. 
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3 Final remarks 



In this work we have analyzed the impact of a general set of soliton solutions over the reflec- 
tion coefficient of the bulk particle collisions with a 3d domain wall, as originally proposed in 
PP. We have shown that when the potential parameters are such that A = 4//, the effective 
potential interacting with the bulk particles have its reflection coefficient arbitrarily larger, 
depending on the value of the degeneracy parameter cq. A remarkable situation is that when 
Co reaches its critical value. In that case, the effective potential becomes a kind of step 
potential, so that the bulk particles having energy lesser than that of the step potential will 
always be reflected, so producing a quadratic driven acceleration. The degenerate structure, 
however, in principle would be created with an arbitrary Cq, far from the original critical value 
of the universe, so leading to an exponentially growing driven acceleration. Therefore, while 
the universe is cooling due to its expansion, the critical value Cq(A(T), /x(T)) approaches c 
and, as a consequence, the reflection coefficient increases and the acceleration also becomes 
quadratic. Thus, supposing that both types of topological structures were created in the 
origin of the universe, one can see that the degenerate structure would dominate the be- 
ginning of the process. This happens due to the fact that, when the reflection coefficient is 
close to one, the scale factor a(t) ~ t 2 and for an exponentially growing reflection coefficient, 
it grows exponentially too. Therefore, in the beginning the degenerate structure dominates 
the expansion and, as the time goes by, both structures predict a quadratic behavior, which 
is compatible with the present experimental bounds. 
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Figure 1: Typical profiles of the soliton solutions in the case where A = \x (a = 1) as a 
function of x, both when c is close to its critical value (c = —2 in this case) (solid lines) 
and far from it (dashed lines). The x field is represented on the thick lines and on the thin 
ones. 




Figure 2: Energy densities for the case where A = 4/x, where cq is near its critical value (solid 
line) ( Co = 1/16 in this situation), and when it is far from this value (dashed line). Here it 
is evident the appearance of the double wall structure. 
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Figure 3: Comparison of the domain wall effective potentials for the case of Brito et al. 
(thick solid line) and degenerate solitons when cq is far from its critical value (thin solid 
line), near it (thick dashed line) and at the critical point (c = —2) (thin dashed line). All 
of them are depicted for the case where A = \i. 




Figure 4: Comparison of the domain wall effective potentials for the case where A = 4/x. The 
thick solid line represents the critical case (co = 1/16). The other curves correspond to Co 
far from its critical value (thick dashed line) and near it (thin solid line). 
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